


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1983 


Development of a microcomputer coupled 
Atmospheric and Oceanic Boundary Layer 
prediction model. 


Tarbet, Gary Lee. 


Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/19748 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
۹ D U DLE Y research materials and institutional publications created by the NPS community. 
dl | Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
KN OX appointed — and published — scholarty author. 


"Wy LIBRARY Dudley Knox Library / Naval Postgraduate School 


411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 








http://www.nps.edu/library 


1٩» با‎ ۰ 
t i* ` eh i 
D 


KM MED? 


4 


ee faa و‎ 


a 
ج“‎ 


"ÉL 


B Se 1 
SRIF) ose 5 


d CNL 0 


5 


AT Lact r, 


EE 


ET: E 
0 یمه‎ 
dr "ke 


` ` 


rated خی‎ + 
did X. 


ن 
BS = 14 d‏ 0 
Cd Le web 4 MET‏ 
CE ile (KIK‏ 


=» m 


Ay بای( وه‎ 
ES "E a SE 
۱ ^ d A APE Ma 
Ni 4 m EE WK 1 AE 
»" me ۸ ۸ 
CA E 5 MEM 
K (TO EET 
e ٩۱ weh 3. A.1 
۷4 a m 


itk PE RR 
2 Ex ez? ا‎ en. 


NICK) Wë EA ET 
ed QU LP ree IST AE 
ae Nar Ag Lyme nn 
۱ ون‎ CALTE P DN 
A ova x la 


H 
ei TI N. +: 
ts ۷۳۴۰۱۵۱ ۸ 
ای رال‎ LT) 
Lr TEE WAN 
A PAR SUD ME 
Joie Kent A 


A CG KEN 
DN NW 
ER 
Aa Ge VER 
GEN 
KE) BE a Ap 
bé Ch KW 41 um T T 
NN. Wa ck 


oo 
E [L4 a 
SH H Je Macs; dl, 
ها ها‎ LEITET 
* De E WM UE? e Kun 


ER‏ نا 
xS ETA‏ 
Ark i‏ 


E 
i LINES 
SAM vr 
Mor de 


Lv 
a? 3 Gs 4 


kW 
Gr 


ck 
Ra KI 
ux Tn 


H Ke Ac ME KO Ee za 
A4 EI LU CHI 
BZ ke تم‎ el 
2 ۱ درمتت‎ 
Lei N LA KAARTEN 


ESA 


IA 
ELS A A $ 
REEL 


ON DON mro 
20 IDE PR ۲ 
AK A KEEN 
s M NR WI 
KH اه‎ 
e ۱2۲ 
>. Siam ss HI 344 4 
9 


۷ مھ UTE‏ "بي با نید 
EA Sark, OREO‏ 
ei ۸ ws a kos ea ee‏ 


LI L'A XE] Kai 1 
Wa Eih) E "Be nd 


9.۰۱ 


H 


Deag 
EP RER SN 
MEE Déi = 


RIA WEI Ka an 


ve b Kë 


KZ Lun 


Fand. Cs." Bal 
3 ۱ E A: 4 
1 e 
d E TM Wu ۰ 

۱ " gi 


A 
bac e A 


۰ 


Ze 3 
E 


WES 
Moni! 


y 
UA PIS 
9 ۱ re 
Rete PIED oe OU PL te 
a. P ی اش در‎ cnr eg Fk Arnd 
MIA E E a عو و‎ a) 1 De وا‎ ۶ 
FM قر‎ E n pl 
es ie 
Pu 
DR 
air, 
GU. ترا ور‎ 
Po 1 
: De gé Pis 
i WS 
WA, UM 1 
cu ra E 


D 6 


u 
4 ye e all Ay 
Nie و‎ pu, 


0 Y ^. Wa 
EE D SE 
dj D Vt wé An 


UL VT 


o d NU 
be LAL 
is MA SH H 3 


H 
Die v Ary d E E E 
u gé fu 
e j E WC Fw 
rh Gabi, 
di Cf DS 


MH 

ko? 

$ po Dn 

i u PAR Gg 7% ër A i 
ربص نا‎ d 

a mrs E 2 2 OI 


n 


BETAN Kar 


j LX MD rah 

NZ ET 

del EN ce AN ie 

A DCH LE Ed ye 

ig 2 SÉ? Mp cm A TAG IN 
LA TUAE pan Un 

E D HP 
A EE: 


rt 


E isa TUM 


P Hue: m 
Y WT È ۳1 


rf 
EL urn fi 


haa 
a LS ADD 
Su St 


^ H Më: Ce 
SE DO ja ia 


ee 3 
Ki; X Lot ery Ee Z 


BLAN ON A D A 
ORT eee Hy y! 
tee LE 


٩. ٩ Ke Wi Je 


SS 


EN, RA 


PM و‎ 


WATER 

Bar‏ ا 
Para a ere Te‏ 

DLP Lei 

A: M E EL Pp 


pants 


Wi np ns "m es 


MTS 
Am EE 


-—. ml, AB 

ELCH te D 

WE 3 

RR KEAA 

Ge A LE c 
ER 


AR Sn x 
Kee um ee 
UD EE 
- ممه‎ 


Se Mott Nue All L4 
LLL G 
TEEN A 

MAIN ۱2: Y 7 

We Ba? Se 


Méi EE $ 
ALLEN i LAU 


ny uvm CP n eit 
ر‎ re Rh BE SET PR 
اد‎ ien LEX 
MS ae 
A. v 
MRT ER ue 
AE Oe مر‎ as CH 
تب‎ rae EE 
eri, EDD h 
n,m Lr EA 
BS) 
NL Aacht MECH A 
EI > ye AT Aoc EXON 
Pic LT born; 0 ۱ Bo, IP 
» A ELIT LCS e Ln u 3 
E D هی سا‎ Caen, Aa Tee 
DL H 


0 
ما‎ 5 = 6١ ch 
SE sn 
۱ USE Se 
۳ ابش‎ bh HA 
KT AE " 2 eh Ai e ^um 
Li je i4 ۰ 4 
MD S ELI 
PO دا‎ Reh A EI 
€ CR E 


Etat 


۳ ORE 
de EIC Dr WA. 
A EN 0 
E 5 EA 
1 ۶ ,مه‎ TM MM erem DN KH le 
b: E AR Dog SE 00 ا و‎ 
iac y E Ki 4 ايو‎ 
ihe hy T چا یه‎ ch SE 
EE 
1۳ ERN ey 
E ER AR DN is Zeg? Ka 
D Ke Re? 
y 


if Sera 
HN Deh E = 
WI 17 
s 3 ee وه‎ 


d 
۷,۹۸۰ و‎ NER df 


E ۰ 


AN Pope ate ۰ 


RER K MEAM 


fn ۳ یط‎ e 


i ak ut ke 


E EES 
NE S ur 


d ENK Kä 

A AMB و‎ 

4 PIU S 
ERK SEG 


a et 


ek GC 


4 


ER Meis Ge 


Se, 
اا ف ا ےک‎ 
ey SE Ne 
wy iw 7 9 
gd E a give an Dune 2 يم‎ LA m e E di Sec 
A En En 2 ره یه‎ ina i € GE? Dr "y d LI 
wi 


sth 


2 
Aide Deet 
^ 


و شرا rA‏ 5 


ES 

V on 

Kë ^‏ یی 
INA, SV Es‏ 

FF qe EQ y 

F II Va 


: bris 
Ae Ee 
Bees d 
GE 


LE aed e 


Ent 0 


EE d Ed keap. 
rd 
"e KEA Lo LE E 
vr hs CN PE 
del, Mey E. pa d MED "e 
بت‎ Ia E d A ir de Der E 
1 DURER Lor E 
Was sh TES e qe 


Ku d dude diu Déi 
Iud SC i. 20 
oe u 
7 er a at 

d‏ ی 

m P ce SE 

E 


E Wé 
dE 2 
"eta 
ATA 
N > iud 


e d E ER 


D 
Es 
Cope 

MS 

30 H ge 
E SE? D dE re 

و 

dte AY a n Es £ EE 
ELA E nen a'y 
rel Cdk epa EA en. cera SE 
D GC A 3 A 
UA yu ee DE 


2 Pd EEE E 
La EA E - Ped A We SEH 
CR AE 


x^ ga d: [t 


dal se gn e i. 


Kéi 0 


A 


E Wa Ze 
d HK. EOM 


ند 


١ 1ه‎ 
I ACE 


PM GEN 
9 den P pend ap Le di 
Ta D E 


KM u 
a i EL 


d d T Pi 
ا‎ O LI H 
EE DEADLY dee Wan 
A A el eme Kay eo 
9 SE un 
kiy E T 
m 
1 Deor NE Ce el, et E 
ES a C4 vw wat 
e Li. 


nahe 


T 


PRA 


VE ET د‎ E 

3۵ y, d ER Pr 

CO? Le sa A 
KR 


ad as eda 
o A Ge 
Lee 
LE 
1 e 


ee 
GE E 
EE 


E 


TE 
p P E. 
RS eee Ze 


e Wu MIA 


0 


ی 
S US‏ 


e Ge 
GE e 
EE 





DUDLEY KNOX LIBRARY 
NAVAL POSTGAADUSTE SONOOL 


MONTEREY, CALIFURNTA “3943 








NAVAL POSTGRADUATE SCHOOL 


Monterey, Galifornia 





THESIS 


DEVELOPMENT OF A MICROCOMPUTER COUPLED ATMOSPHERIC 
AND OCEANIC BOUNDARY LAYER PREDICTION MODEL 
by 


Gary Lee Tarbet 


December, 1985 


Tesis Advisor: | Kenneth Davidson 


Approved for public release, distribution unlimited 


1218054 





Val LA 


DUDLEY KNOX LTPPADY 





SECURITY CLASSIFICATION OF THIS PAGE (When Dete Entered) NI 


REPORT DOCUMENTATION PAGE e de ی‎ 


. REPORT NUMBER 2. GOVT ACCESSION NO. 3. RECIPIENT'S CATALOG NUMBER 

















4. TITLE (and Subtftie) 5. TYPE OF REPORT & PERIOO COVERED 


Development of a Microcomputer Coupled 
Atmospheric and Oceanic Boundary Layer 
Prediction Model 





Master's Thesis 


December 1983 
6. PERFORMING ORG. REPORT NUMBER 











7. AUTHOR(a) 8. CONTRACT OR GRANT NUMBER(S) 












Gary Lee Tarbet 





10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 





9. PERFORMING ORGANIZATION NAME AND ADDRESS 






Naval Postgraduate School 
Monterey, California 93943 










1. 12. REPORT OATE 







CONTROLLING OFFICE NAME AND ADDRESS 








Naval Postegraduate School „Decenbe 98 
Es 2 ” eg ۰ FPA 
Memcerey, California 93943 ۳ = 


- MONITORING AGENCY MAME & ADDRESS(i1I different from Controliing Office) 15. SECURITY CLASS. (of this report) 


iSe. OECLASSIFICATION/ DOWNGRADING 
SCHEOULE 


۱6۰ DISTRIBUTION STATEMENT (of thía Report) 


Approved for public release, distribution unlimited 





17. DISTRIBUTION STATEMENT (of the ebatract entered In Block 20, lí different from Report) 


| 18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reverse aide if neceseary and identity by biock number) 
Marine Atmospheric Boundary Layer, MABL, Oceanic Boundary Layer, OBL 


20. ABSTRACT (Continue on reveree side if neceseery and identify by block number) . 
A coupled Marine Atmospheric Boundary Layer (MABL) and Oceanic Boundary Layer 


(OBL) model is developed using the Naval Postgraduate School and Garwood 

models respectively. All coding is done on the Hewlett-Packard 9845 micro- 
computer with emphasis on ease of use. The model is used to explore cases 

when feedback between the boundary lavers significantly influences model fore- 


Ecc ie Sensitivity of the model to slight input variations is explored. | 
Light wind situations where stratus or fog formation is extremely difficult to! 
predict is investigated. Cases covered include variations in mixed layer depth 
and wind speed which produces significantly different forecasts from the 


2 TEM 


DD , ee 1473 EDITION OF ! NOV 65 IS OBSOLETE 


S/N 0102- LF- 014- 6601 SECURITY CLASSIFICATION OF THIS PAGE (When Data Entera 





Approved for public release, distribution unlimited 


Development of a Microcomputer Coupled Atmospheric and 
Oceanic Boundary Layer Prediction Model 


by 


SAS Se 
Lieutenant Commander, United States Navy 
Ens University of Utah, 1975 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN METEOROLOGY AND OCEANOGRAPHY 


from the 


NAVAL POSTGRADUATE SCHOOL 
December, 1983 





ABSTRACT 


A coupled Marine Atmospheric Boundary Layer (MABL) and 
Oceanic Boundary Layer (OBL) model is developed using the 
Naval Postgraduate School and Garwood models respectively. 
All coding is done on the Hewlett-Packard 9845 microcomputer 
with emphasis on ease of use. The model is used to explore 
cases when feedback between the boundary layers significantly 
influences model forecasts. The sensitivity of the model to 
Slight input variations is explored. Light wind situations 
where stratus or fog formation is extremely difficult to 
predict is investigated. Cases covered include variations in 
mixed layer depth and wind speed which produces significantly 


different forecasts from the initial input. 





E. IDE ND ALO a. EN. ... + 5 sls De TN. vs 

TI. DESCRIPTION OF BOUNDARY LAYER FEATURES 2 ۰ . ۰ ۰ ۰ » + 

۱۳۳ MODELS ...... Joaca o e... e esee se فلل بك‎ 
A. MARINE ATMOSPHERIC BOUNDARY LAYER (MABL) 

A "LUE UNUS... see. 
ENEOCEANUECOBOUNDARY LAYER (OBL) MODEL ....... t 
LCSNUESSDISOUMDARYSEAYERCMODELDL. .....e999992992* 

IV. PN NUM EESULTS ELS. Que ema e" eee ف‎ 
BRETT ESBERED MODELFZRESUETS, 2.2.2.2... e... ...... 
e OIDO EA ERA 
C. COUPLED MODEL VARYING MIXED LAYER DEPTH 

E LE UNE. e. s.s» on 

Pre DOUBLED MODEL VATYING WIND SPEED RESULTS ...... 
V. RESUS AND CONCLUSIONS در و ههد وه‎ SEN 
net ON NEN LS هو‎ ۰ ۰ a عع .مم‎ SOT 
II eege ee هه‎ o e sn e e 


TABLE OF CONTENTS 


16 


16 
23 
28 
355 
35 
D 


41 
44 
48 
50 


32 





IU. 


lii"; 


NP. 


ES 


14. 


LIST OF FIGURES 


Simplified Atmospheric and Oceanic Boundary Layer 
Temperature Profiles e © © ه©‎ 0 ée e © © © © © © e © e © © © © © © © © © © © © © © e 


Mechanical Energy Budget for the Ocean Mixed Layer 


Simplified Flow Diagram of the Boundary Layer Model 
Showing Possible Configurations of Input Informa- 

tion, Interrelation Between Atmospheric and Oceanic 
Models, Model Outputs, and Tactical Models Which Use 
Moe TOME DWES .....99 ne IMMUNE. Sr S s 


Input and Flow Chart for MABL Prediction Model ..... 
MEBE and Elow Chart for OBL Prediction Model ...... 


Input and Flow Chart for Coupled OBL and MABL 
Prediction Model ...... ی‎ o's. is oes w vss 


Data Set e e e e e e © © 0 e e © © © © © © © © e © 0 è © © © © © © © e © © e e e e © © © © © © 


Uncoupled Air Sea Boundary Layer Model 24-Hour 
85152095256 Using Fixed SST and MLD Values ............ 


OBL 24-Hour Forecast for Original Input Conditions 


Coupled Air Sea Boundary Layer 24-Hour Forecast 
ADO ECON EOS um wm ew e el NNN 


OBL Forecast With Initial MLD Set at 10 Meters 
AMOCO Aput Values are Held Constant ........... 


Coupled Air Sea Boundary Layer 24-Hour Forecast 
Varying Initial Sea Surface MLD to 10 Meters Vice 
۱۳۳۰۹ LEO 2 ۱16۳6۲ و مد و و ود ت‎ nenne 


OBL 24-Hour Forecast With Wind Speed Increasing 
CORO Es” In the Forecast Period ...:....-.-.- 


Coupled Air Sea Boundary Layer 24-Hour Forecast 
Varying Wind Speed from 4 to 10 Knots in the 
ORECAST Perlod Menn. Soo 


11 


Le 


14 
18 


29 


351 
34 


36 
38 


39 


42 


43 


45 


47 





I. INTRODUCTION 


Military leaders, engineers, and scientists have become 
aware of the environmental effects on electromagnetic (EM) 
emomelectrooptical (EO) signal propogation. Many of the 
current weapons guidance systems, command and control 
communications, and electronic countermeasures are critically 
dependent upon environmental parameters. One extreme case 
was recently brought to the attention of Pentagon officials 
when a new missile guidance system became totally ineffective 
in certain environmental conditions. The modern naval leader 
must not only be aware of the environment but must also know 
how to use the current environmental conditions to best 
advantage. The deployment of resources, decision of 
appropriate weapons systems, and overall tactics must include 
a consideration of EM/EO propagation. The overall 
effectiveness and the successful outcome of an operation 
could be tied to this very knowledge. 

Dnncesctmospheric factors which effect EM/EO propogation 
are the temperature, humidity, vertical gradient of pressure, 
small scale inhomogeneities or turbulence, distribution of 
aerosols, and concentration of water vapor. The refraction 
of EM/EO signals is primarily affected by the first three 
factors.  Turbulence affects the index of refraction through 


wave front distortions while the remaining factors cause 





Ebrunctaron and dispersion. All of these effects are 
interlinked and they must be computed simultaneously. 
Another problem comes in measuring each of these factors. No 
known or planned system provides the accuracy required for 
direct measurement. Indirect methods will have to be employed 
by units in the operational arena for the foreseeable future. 

An equally important problem is prediction of the mixed 
layer depth and sea surface temperature in the ocean. 
Research is underway to find the relationship between the 
synoptic scale weather patterns and the sea surface 
temperature. Naval operations generally take place in areas 
where the marine atmospheric boundary layer (MABL) has been 
extensively modified by contact with the ocean surface. The 
effects of heat flux from the ocean to the atmosphere warming 
the boundary layer, the subsequent increase in turbulence, 
the transfer of water vapor to the air, and the effects of 
Salt and other aerosols being injected into the air by waves 
are all important to Naval operations. These fluxes of heat 
and water vapor can change the structure of the MABL to the 
ExLent that clouds or fog are formed. Clouds and fog will 
dramatically reduce the short wave solar radiation striking 
the ocean surface thereby reducing the surface heating due to 
radiation. The diurnal cnange in the sea surface temperature 


and mixed layer depth will be decreased. 





In determining the effect upon acoustic propagation, both 
the depth and strength of the mixed layer gradient must be 
considered. Skip zones and ducting are examples of oceanic 
phenomena which must be considered in every naval operation. 
Unusually strong or weak diurnal affects can significantly 
alter these factors. To provide an optimum forecast of the 
OBL, the effective shortwave radiation, internal mixing 
forces, and atmospheric entrainment must be taken into 
consideration. 

It is obvious that any attempt at modeling the MABL and 
OBL should be linked for optimum results. Microcomputer 
programs developed at the Naval Postgraduate School for the 
MABL (Davidson, et. al.) and for the OBL (Garwood) have been 
linked to provide the necessary feedback. While both of 
these models have been verified independently, the linking 
should improve forecast accuracy. 

Determining situations where the linking has significant 
effects rs the primary goal of this thesis. 5 22515121 7 
studies were also conducted in an attempt to determine which 
if any of the factors provide significant differences between 
the coupled and uncoupled models. Additionally, since the 
Model is currently running on a Hewlett Packard 9836 
microcomputer and the fleet units are and will be using 
Hewlett Packard 9845's for several more years, it is 


necessary to transfer tne code to the latter unit. Since the 





internal architecture of the two systems is significantly 
different, changes in program structure will have to be 
verified for accuracy and consistency with the original 
model. 

Having mixed layer forecasting capabilities onboard 
should enable the operational fleet units to use the 
environment to maximum advantage. Not only can forecasts be 
updated rapidly as on-site conditions vary, but those "what 


if" questions can be answered quickly and accurately. 





II. DESCRIPTION OF BOUNDARY LAYER FEATURES 


The MABL extends from the surface through the capping 
inversion which is typically .5-1.5 km above the surface. 
The MABL is cooler and.more moist than the overlying air, and 
it is capped by an inversion 50-100 meters thick. 
Temperature increases and humidity decreases with height in 
this inversion. The air-sea interface is bordered by oceanic 
and atmospheric turbulent mixed layers which effectively 
M late the quasi-geostrophic regions above the inversion 
and below the thermocline. The OBL or mixed layer in the 
Ocean typically spans the upper 10-100M of the ocean. Mean 
velocity and density values tend to be vertically uniform in 
this region. At the bottom of the mixed layer a transition 
region exists called the thermocline. Turbulence in these 
well mixed regions is created by bouyancy, flux and velocity 
gradients that are a result of air-sea interactions. The 
vertical homogeneity of these two mixed layers can be 
attributed to the strong mixing by the turbulent motion. 

Bouyancy driven energetic eddies fill the OBL and MABL. 
In the atmosphere the eddies entrain warm, dry air with high 
momentum from the free atmosphere into the boundary layer. 
If this entrainment causes the MABL to extend above the 
lifting condensation level, then clouds or fog will form. A 


typical profile of the MABL and OBL is shown in Figure l. As 


10 
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Figure 1. Simplified Atmospheric and Oceanic Boundary Layer 
Temperature Profiles 


abt 





ET De seen from Figure 2, the bouyancy driven fluctuations 
have an even more direct role in the mechanical energy budget 
for the OBL. 

With the understanding of the importance of bouyantly 
driven entrainment effects, the necessity to couple the near 
surface prediction models is ebene A cause and effect 
relationship is developed through the interactions of the 
ocean and atmospheric surfaces. Examples of this effect 
include: 

sur fscerbouyaney flux induced entrainment not only 
increases the depth of the MABL, but it also changes its 
effect on the ocean mixed layer. 

2. Clouds in the MABL can be caused by a change in the 
ocean surface temperature which in turn affects the 
radiation budget. 

The relatively complex models used to predict these 
features have been tested in both the coupled and uncoupled 
modes. While often little improvement is noted in the 
output in the coupled versus uncoupled modes, under certain 
circumstances the coupled mode is mandatory and produces 
"0 0111121232197 better forecasts. It is the goal of current 
research to determine exactly what factors have the strongest 
influence on coupling and under what circumstances coupled 
models must be utilized. 

The purpose of this thesis is to show under what 


circumstances the coupled approach is most useful. The 


answers will hopefully be obtained through interpretive 
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Figure 2. Mechanical Energy Budget for the Ocean Mixed Layer 


efforts uSing the power of the computer. By varying the 
angle of radiation (latitude), amount of radiation (month), 
and start time for the program, each set of output will be 


compared and analyzed for consistency. A typical case will 


then be selected and further studies will be conducted. The 


effects of coupled versus uncoupled oceanic mixed layer depth 


Variations and wind speed variations will be examined. 
the magnitude and usefulness of this effort is 


illustrated in Figure 3. Outputs from the coupled model can 
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be used tactically in forecasts shown in the righthand column 
of Figure 3. The models to-be used for the forecast include 
the Garwood Model and the Naval Postgraduate School Marine 
Atmospheric Boundary Layer Model. A detailed description of 
each model and the approach taken in coupling the outputs 
will be discussed. While this understanding is not necessary 
Bene Utlltzation of the output, it may provide useful 


information in obtaining maximum benefit from the program. 


3 





III. MODELS 


A. MARINE ATMOSPHERIC BOUNDARY LAYER (MABL) MODEL 

The NPS MABL model is a zero-order, two layer, integrated 
mixed layer model. The model assumes the atmosphere consists 
of two layers; a well mixed, turbulent boundary layer, and 
the relatively non-turbulent free atmosphere above. The 
model is based on radiative transfers described by Davidson, 
et. al. (1983) and entrainment energetics formulated by Stage 
eae Businger (1981). 

The two zones are separated by an inversion layer or 
Eransition zone. The zero-order model assumes this 
ans ition zone to be infinitely thin; therefore, a jump or 
EN-continulty occurs at this point in the profiles of all 
conservative parameters. 

The current model requires the following inputs: 

l. An initial atmospheric sounding. 
2. The geostrophic wind. 
3. The surface temperature. 

In an operational scenario the boundary layer winds can 
be estimated from standard meteorological charts. During 
this evaluation the actual hourly winds were input for 
initialization purposes. In the uncoupled version, sea 
Surface temperature (SST) remains unchanged; however, when 


coupled, the SST does change with time and is predicted by 
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the OBL model. The humidity and temperature in the well 
mixed MABL are predicted. Inputs of the surface wind and 
wind shear at the inversion are required by the MABL model. 
The large scale subsidence normally obtained from synoptic 
scale NWP products must also be prescribed for the model 
period. 

The model is very sensitive to subsidence values; 
therefore, care should be taken in selecting this value for 
proper results. Three methods can be used to compute the 
subsidence (large scale vertical velocity) from single 
station observations. These three methods are the kinematic 
method, adiabatic method, and integration of the moisture 
budget equation (Q-method) which are all well described by 
Gleason (1982). Gleason's study showed the Q-method 
displayed the most merit as a single-station assessment of 
subsidence. Computation of the solar zenith angle, which is 
used to compute effective short wave flux, uses the latitude, 
julian day and start time which are initial input values. 

As shown in Figure 4 the atmospheric model has a 30- 
minute time step. During each cycle, the program predicts 
the mixed layer temperature, humidity, and the jump of these 
values at the inversion. When clouds or fog are formed, the 
cloud top cooling and entrainment computations are important 


in the physical processes. 
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Prescribed inputs: 
Aeostrophic wind speed 
sea Surface temp. 
Subsidence 
Advection 


le rofiles: 

emperature 
Humidity 
(Surface to 3Km) 


te: 
surface fluxes 
Condensation levels 1 





Compute: 
Entrainment ntrainment  . 
Croud "top cooling 


Compute updated: 
Mixed var depth 


Mixed layer temperature 
Mixed layer humidity 
Jump strengths 











forecast 


complete 
Output: 
Graphs 


Figure 4. Input and Flow Chart for MABL Prediction Model 


EM ducts 


Optical turbulence and extinction 
Dispersion 
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The conservative quantities and their jump at the 
inversion is predicted by using the standard integrated rate 


equations (Tennekes, et.al., 1981). The equations are: 


h(Dx/Dt) - zn - Ga + source (1) 
h(Dax/Dt) = hrx(ah/at) - (w'"x")pg + (w'"x")h ¬ source (2) 
۳ = Lapse Rate 

MU rec zz tene 

m * Net Radiative Heat Flux 


۱۱ subscripts "h" and "O" refer to inversion height and 
surface values respectively. 

To close this system of equations and to compute the 
variation in the inversion height (Stage, et.al., 1981) 
entrainment velocity parameterization is used. One 
additional assumption is used to close the system and that is 
that the dissipation rate of turbulent kinetic energy (TKE) 
is a fraction (1-A) of the production rate. The entrainment 
coefficient (A) is taken as .2 for the formulation. 

The bulk aerodynamic formulas are used for surface fluxes 


of momentum, sensible heat, and latent heat. 


ECE T (3a) 
m= e ilog - o) (3b) 
SES Gerd ten - q) (3c) 


These fluxes are given by: 
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u'w' = u*2 ` (momentum) (4a) 








T'w' = u*T* (sensible heat) (4b) 

q'w' = u*q* (latent heat) (ie) 

ey and - ten meter stability dependent drag 
coefficient 

¢ = potential temperature 


Speceitie humidity 


iQ 
1 


The subscript O0 denotes surface values. 

In the MABL model, extensive work has been done on the 
۰ ۱ 1 1 0۳۱ portion because of its significance to the OBL 
model and the sampling of the two models. The short wave 
radiative flux is computed through use of the delta-Edington 
Method. An excellent review of the delta-Edington Method 
including all parameters, atmospheric factors, and equations 
has been published by Fairall (1981). This portion of the 
model was added to account for the heating of the mixed layer 
by solar radiation. 

In the boundary layer, short wave extinction is dominated 
by scattering vice absorption. This second short wave 
radiative component is usually referred to as diffuse solar 
radiation. Atmospheric particles such as cloud droplets and 
sea~salt aerosols are the primary scattering nuclei in the 
MABL. The current MABL model computes both direct and 
diffuse radiation components to determine a total short wave 


radiation flux value at the surface. EPnsasdadıtion, suche 
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fraction of reflected short wave radiation, Ag, from the sea 
surface is prescribed as .l in the MABL model. 

The modeling of radiative flux transfer has been 
accomplished in numerous ways; however, even the simple 
models are extremely complex compared to other 
parameterizations used in the model. Even though long and 
short wave radiative fluxes are computed Separately, there 
are numerous sources of error in these calculations. Some of 
these sources include: 


WwW Concentrations Of lack of absorbing gases such as 
carbon dioxide, ozone, or water vapor in the 
atmosphere. 


EX Theouncertainty of quantity, size and distribution of 
background aerosols. 


EE hec Size of various cloud droplets and their 
distribution. 


Since this model primarily intended for use over ocean 
areas non-black stratus clouds were permitted by introducing 
mul emissivity (e.) into the long wave radiative flux 
۳۳۲ ۱ ۱] 1 017۰0 Cloud emissivity is a function of total cloud 
۱۳/۳۲ ۱0 content, w. Cloud liquid content profiles are 
Ecouoxrmatelyv linear with height (Davidson, et.al., 1983) 


Cloud water content and emissivity are given by equations 5 


and 6. 
w = 0.5 p. (h-Zec) Gh ey) 
a = l - exp(-aw) (6) 
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alr density (1.25 X 1073 gm cm^?) 


Pa 
h = height of the mixed layer (cloud top) 
Pr being eondensatron level (cloud bottom) 


BE ado, et. al... 1991) 


Gp, = liquid water content at cloud top 


Using the Stefan Boltzman Law, the net long wave cloud 
Eop radiation flux, Lnpp» can be calculated from the cloud top 
temperature, Tp. The cloud bottom temperature, Tc, and the 
Sea surface temperature, Te, are used to calculate the flux, 


iat the bottom of the cloud. These fluxes are given by: 


Inh = eco ) = EIERE (7) 


Mur To) (8) 


1 


L 


E 


nc C 


Stefan's Constant (4.6l x 107i 


Q 
1 


Eq = obtained from equation (4) 
The net long wave radiation at the surface, Flong' becomes: 


MENT. - ۰. = (1 > ec) Tay) (9) 


T = average cloud temperature 
For the cloud free case, the net fluxes are calculated at 
ana O0 by integrating the flux emissivity profile 


(Fleagle, et.al., 1978). Ine 5 ٩ ۵۲۱8 wave flux at the 


Surface for the clear sky case 1S given by: 
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Fiong = Fu - Fa (10) 


UDI radiative flux‏ أن 


downward radiative flux 


Fa 


B. OCEANIC BOUNDARY LAYER (OBL) MODEL 

A mixed layer model for the ocean using the continuity 
equation for an incompressible fluid, the first law of 
thermodynamics (heat equation), the conservation of salt 
equation, the Navier-Stokes equation of motion with the 
geostrophic component eliminated, an analytical equation of 
State, and a two-component vertically integrated turbulent 
kinetic energy budget was developed by Garwood (1977). 

An understanding of the dynamics of the entrainment 
process is a key factor in predicting the variable changes in 
the mixed layer. The turbulence of the overlying mixed layer 
provides the energy needed to destabilize and erode the 
underlying stable water mass (Garwood, 1977). The turbulent 
Kinetic energy equation is the basis for the entrainment. A 
closed system of equations is obtained by using the bulk 
buoyancy and momentum equations with the mean turbulent field 
modeling of the vertically integrated equations for the 
individual turbulent kinetic energy (TKE) components. 

To better define the mixing process, separate horizontal 
and vertical TKE equations are used. Energy for vertical 


۲۳۱ تا‎ IS provided by both buoyancy flux and shear 
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production. The buoyancy equation is derived from the heat 
and salt equations coupled with the equation of state as 


shown in equation (11). 

og [1 - a(8- eg) * &(S - so)] GE‏ = م 
Beuyancy is given by:‏ 
b= g (0g - (0 (12)‏ 


9 = temperature 
s = salinity 
p = density 
g = gravity 
a = expansion coefficient for heat 
g = density coefficient for salt 
Note: The tilde represents instantaneous values and the 
subscript 0 represents an arbitrary, but representa- 
tive, constant value. 
The effect of the salinity on the short-term density profile 
evaluation is generally found to be insignificant except at 
higher latitudes. Temperature is usually the dominating 
factor in the density profile. However, by using buoyancy 
instead of only temperature permits the model to be applied 
in Situations where evaporation and precipitation contribute 
Significantly to the surface bouyancy flux. 
For extended forecasts, the Ekman wind-driven horizontal 


current profiles as well as the temperature and salinity 
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purses must be provided with initial values. The mixed 
layer depth, "h", as defined by the Garwood OBL model, is the 
shallowest depth at which the 1 density value, SEES 
۲۱۲۲ Units greater than the observed surface density value. 
Additional ocean parameters which must be prescribed include 
Buemradration extinction coefficient, the fraction of short 
wave radiation absorbed in the upper meter of the ocean, and 
Mmmemeritical Richardson number which defines a stability 
adjustment at the bottom of the mixed layer. Surface 
boundary conditions required for the OBL model include air 
temperature (dry bulb), dew point temperature, wind speed and 
direction, the rate of evaporation (E) and precipitation (P), 
۶۳۲ 1 incident solar radiation. 

Using the bulk aerodynamics formulas, the turbulent 
les of sensible heat, Qn, and latent heat, Qe; Can be 


computed as follows: 


Os = Cg (.98 Eg - E4)U10 (133) 


Op = Ca (Ts - Ta) Vio (13b) 


The net back radiation is estimated from the empirical 


dauat ron (Husby, 1978). 


Op-1.14x10 7(273.16+7,,2(.39-.52,1/2) (1-.6c2) (13c) 
2 Saturated vapor pressure (.98 corrects for salt 
defects) 


Ba = vapor pressure Of air based on dew point temperature 


نات ۳ کل ۱۳ 
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sea surface temperature 


H 
u 
il 


fractional cloud cover 


Q 
I 


Geewupward heat flux, Qu; 15 then given by: 

20 = Oe +t Oh + Op (14) 
The solar radiation, Qs, iS given by: 

(15) و (6663, - 1( P)‏ ج - 1 < ۲ 


The constants "a" and "b" are adopted from Tabata (1964) and 
the cubic cloud cover correction from Laevestu (19600). 0015 


the clear sky radiation given by Seckel and Beauday (1973): 
مث - إن‎ ۲ ۸۱ coso + B] sing * A9 cos29$ + و85‎ 5182 4 (16) 


EN ostficilents AQ, Aj; etc. are calculated by harmonic 
representation of the values predicted in the Smithsonian 


Meteorological Tables with 
SOS) (t—21) (17) 


where t is the julian day of the year (O'Loughlin, 1982). 

A very small percentage of the incoming solar radiation 
penetrates the ocean mixed layer. Approximately 50 percent 
is absorbed in the first meter of the ocean in most parts of 
the open ocean. The portion absorbed varies from region to 


region and is highly dependent upon such things as suspended 
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particulate matter and phytoplankton. More radiation will be 
absorbed in coastal regions than the open ocean because of 
the increased amount of suspended particulates. This portion 
of the absorbed radiation is considered to be part of the 
upward heat flux because very little of this heat is 
entrained into the deep ocean. Most of this energy is 
transferred upward out of the ocean and back into the 
atmosphere. The remainder of the short wave radiation does 
penetrate the mixed layer; however, an exponential 
attenuation does take place which is highly dependent upon 
uter turbidity. WIE The TTacltlOoh of “Solar radiation 
absorbed in the first meter, RF, the net heat at the surface 
E Given by: 

Onet = Qu + (RF) Qs - Qs (18) 

From the equations discussed above, the momentum and 
Surface fluxes of buoyancy (heat and salt) can be computed. 
The mixed layer temperature, Salinity, bouyancy and velocity 


Fluxes are given by: 


۱03 = Onet Cp (19a) 
(S'w'") = (P - E) Sọ (19b) 
A E gie (Tw) = g(s™w")] Cis) 
(u'w!) = U*2 (19d) 


werpt O refers to surface value. The friction velocity 


ai Usas given by: 
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IAS a (20) 
where ts = Pa Cg U 0? (21) 
er = surface stress (dynes cm”2) 
Peel ve Surface buoyancy flux results when Qnet < 0 and 
E > P. During daytime periods, the solar heating at the sur- 
face dominates giving a negative buoyancy flux. At night the 
combination of long wave radiation and the upward turbulent 
fluxes of heat and moisture produce a positive buoyancy flux. 
The ocean model, as shown in Figure 5, details the inputs 
discussed above. At each one-hour interval, new mixed layer 


depth, temperature, salinity, and wind-driven current 


profiles are predicted. 


C. COUPLED BOUNDARY LAYER MODEL 

The advantage of linking the two models described in the 
previous two sections is obvious when examining the inputs to 
each of the models. Allowing feedback of current input 
parameters to occur between the models at each time step can 
Bertenciaitly produce significantly better forecasts. 
5 0010111111 (1982) accomplished the initial coupling on a 
Hewiett-Packard 9836 microcomputer taking care not to alter 
the physical process in each of the models and insuring all 
variable units were passed uniformly. The MABL model only 
requires the SST from the OBL model. The correct SST IS 
extremely important to the MABL model and affects the entire 


output package as discussed in the next section. 
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Prescribed pe: 
Surface wind speed and direction 
Incident Radiation 
Precipitation - Evaporation 
Surface (lm) Radiation Absorption 
Critical Richardson Number 


ired Profiles: 
emperature 
Salinity 
Velocity 


te: 
surface Fluxes 
Stability Parameters | 


l hour 










Step 
Rate of 
Mixing 
function 
(stabilit 
parameter 


Entrainment Retreat 
Mode Mode 


Diffusion 















Updated:‏ زاس 

Mixed layer depth 24-hour Yes 
Well mixed temperature (SST) forecast 
Well mixed salinity complete 





E 
Grapns 


Figure 5. Input and Flow Chart for OBL Prediction Model 
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Niort ele lmcoupling Problems overcome by O'Loughlin 
eT. 13] included: 

1. The atmospheric model uses a 30-minute timestep 
while the ocean model uses a l-hour timestep. The coupled 
model calls the ocean model on every other timestep to 
overcome this problem. 

2. The atmospheric model requires only wind speed 
and not direction. The ocean model requires wind direction 
Compute thewhorizontal ocean turbulent velocity flux, 
U,*^, for the momentum budget equation. A Subroutine was 
added to compute the horizontal wind components from speed 
ma direction input during the initialization. 

A complete flow diagram of the steps in the coupled model's 
prediction computation is shown in Figure 6. 

In 1982 and 1983 the Naval Oceanography Command purchased 
and distributed Hewlett-Packard 9845 microcomputers to all 
aviation support ships and selected detachments. These units 
were designated as interim TESS (Tactical. Environmental 
Support System) units until the TESS system is deployed. The 
9845 has proven itself as a structurally strong computer. 
Many application packages have been written for the unit and 
more are being distributed by The Naval Environmental 
Prediction Research Facility (NEPER)F all the time. 

The original formulation of a coupled model was done on a 
Hewlett-Packard (HP) 9836 computer. Transferring this 
working code from the 9836 to the 9845 would, on the surface, 
appear to be a trivial matter. On the contrary, the 9836 is 


a 16 bit computer system based on the Motorola 68000 micro- 


processor capable of addressing one megabyte of memory. 
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eostrophic 
Solent tail) 
Se 
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Condensation levels 





No louds Yes 


temp 
Mixed layer depth 
MARL, 


Entrainment 
Cloud top cooling 


Campute updated: 
Mixed layer depth 
Mixed layer temperature 
Mixed layer humidity 
Jump strengths 
EM ducts 
Optical turbulence and extinction 
Dispersion 





Figure 6. Input and Flow Chart for Coupled OBL and MABL Pre- 
diction Model 
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Having 16 bit accuracy and the large memory addressing 
capability allowed the relatively easy coupling of the 
initial model. The 9845 is advertised by HP to be a 16 bit 
computer with a proprietary processor to HP. The processor 
has the limited capabilities of an 8 bit processor in 
addressable memory (64K). This required extensive changes in 
the structure of the coupled model. In addition, while the 
company claims 16 bit accuracy with the 9845, the extensive 
changes in code require verification that model physics and 
output have not been modified by the lack of precision or 
round off error within the computer system. Another factor 
in preparing the program was to reduce its overall size so as 
to use only one tape for the program and one tape for data to 
eliminate the confusing practice of continually swapping 
tapes during program execution. Transfer of information from 
Or to tape units is slow and must be limited. Making the 
program as user friendly as possible is an additional 


consideration. 
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IV. DATA AND MODEL RESULTS 


The data set used for this analysis of the model is a 
modified set from the Cooperative Experiment on West Coast 
Oceanography and Meteorology (CEWCOM-76) shown in Figure 7. 
The data set was modified to maximize the effects on the 
program output. Care has been taken to ensure the input is 
reasonable and representative of the area to be discussed. 
The primary goal of this application is to show the 
sensitivity of the program to variations in input. 
Additionally, these model results could easily form a 
Scenario for a fleet application showing the utility of the 
program. 

The problem to be posed for this analysis is "will clouds 
form within the next 24 hour period?" This problem could be 
quite significant if perhaps the forecaster was on a vessel 
with only the local observations. The availability of good 
facsimile and satellite products have greatly reduced the 
burden on present day forecasters. The other problem 
associated with forecasting for an afloat unit, especially 
any U.S. Navy ship, is that these forecast officers transfer 
positions. The forecaster does not have the opportunity to 
gain the expertise of an individual permanently assigned to 
@me forecast office. Therefore, it is of paramount 


importance that adequate tools be made available to the 
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ATMOSPHERIC DATA SET 


Date 21 June 
Latitude رد‎ EN 
Surface Temperature Omen 
Temperature Jump at Inversion SE 
Lifting Condensation Level 567m 
Inversion Level 607m 


يم 


Winds Average 23.5 knots 
Mixed Layer Specific Humidity OZ OG 


Jump Strength -2.4 g/kg 


OCEAN DATA SET 


Sea Surface Temperature 2 ۵۳ C 
Mixed Layer Depth (initial) 20 2 
Jump Strength 2.0 6C 


Figure 7. Data Set 


۳۲ ۳ ۳۳۰ This program is just such a tool. Using onboard 
HP9845 assets, the fleet geophysics officer can input local 
observations and receive a 24-hour forecast for the OBL and 
MABL. In addition this program allows the forecaster to 
answer those nagging and sometimes critical "what if" 


questions such as: 
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l. What if the wind speed varies? 
2. What if the mixed layer depth changes? 
3. What if the subsidence rate varies? 

To examine the coupled and uncoupled models and their 
interactions, the models were initialized using the following 
Semereions. The overall synoptic situation is very stable. 
A large high pressure system is dominating the synoptic 
pattern in the region. Light winds averaging approximately 3 
Knots with the strong subsidence of the high pressure system 
has resulted in clear summer days. No change in the general 
synoptic pattern is forecast for the region by the numerical 
weather prediction (NWP) products produced from Fleet 
Numerical Oceanography Center (FNOC). The air temperature 
has remained around 19° C. The long periods of sunlight have 
caused a strong, shallow mixed layer to form on the surface 
of the ocean approximately 2.0 meters deep with a 29 C 
temperature jump at the boundary. Two soundings and hourly 
meteorological observations have been taken in the last 24 
hours and are available for use. The date is June 21st and 
all data assumes a latitude of 30° N. Start time for each 


forecast is 1900. 


A. Or OUPLCED MODEL RESULTS 
۱۱۳ ۲ اه‎ output to be examined is that of an uncoupled 
(atmospheric model only) model. As seen in Figure 8, the 


mixed layer depth (MLD) is held fixed at the original value. 
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Figure 8. Uncoupled Air Sea Boundary Layer Model 24-Hour 
Forecast Using Fixed SST and MLD Values 
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Winds for this run were forecast to remain light and variable 
and generally out of the north. The specific humidity shows 
an almost linear rise throughout the forecast period as heat 
and moisture is transferred from the ocean surface into the 
MABL. Evidence of this heating can also be found in the plot 


for the MABL temperature. The plot shows an almost linear 


increase toward the sea surface temperature (SST) during the 
first 20 hours. Looking closely, a slight steepening of the 
went does occur with the rising of the sun, and the 
gradient decreases sharply late in the period as the solar 
altitude decreases and the strong temperature difference has 
been removed. The difference in height of the lifting 
condensation level (LCL) and the inversion decreases early in 
the period and then becomes parallel. Looking: at this 
forecast, no clouds will form, however, the LCL and inversion 
height are very close together and asking a couple of the 


"what if" questions listed above would seem appropriate. 


Be COUPLED MODEL RESULTS 

In Figures 9 and 10 the differences in the coupled model 
output can be examined. As discussed earlier, many cause and 
effect relationships exist between the ocean and atmosphere. 
The changes in SST and mixed layer depth are input at each 
time step into the atmospheric model. Changes in air tempe- 


١ 0 ده‎ ana winds are fed back to the ocean model. It would 
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Figure 9. OBL 24-Hour Forecast for Original Input Conditions 
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Figure 10. Coupled Air Sea Boundary Layer 24-Hour Forecast 
EI Initial Input Conditions 


39 





Seem logical that the coupled model should produce 500 
forecast since many of these factors are taken into account. 
The strong gradient near the surface is evident in Figure 
9. Note that a small decrease in temperature at the surface 
has a marked effect on the depth of the mixed layer. Through- 
out the early period heat is being transferred into the 
atmosphere with no replenishment. This trend is reversed 
later in the day when short wave radiation absorbed by the 
sea surface is converted into heat, re-establishing the 
Shallow mixed layer. Changes in the mixed layer are also 
traced in the lower plot of Figure 10. The wind speed has 
been prescribed and is the same as for the previous case. 
The specific humidity curve is markedly different. While the 
early results show the same increase approximately 14 hours 
into the forecast, a strong decrease is noted in the specific 
۳۲۲۲۱ 1 11:۲۰ This change is associated with the formation of 
clouds in the MABL and reduction of the moisture flux as the 
atmosphere becomes warmer than the SST. The temperature pro- 
files are also markedly different. Allowing the SST to vary 
at each time step allows the atmosphere and sea surface tem- 
peratures to come together very rapidly. The sea surface 
continues to cool until short wave radiation inputs reverse 
the trend. Early in the period the LCL and inversion heights 
are Similar to the uncoupled case; however, the feedback 


process does allow the two levels to intersect, predicting 
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EEN-ormatron of a thin cloud deck. In time this cools the 
atmosphere by cloud top radiation, allowing the stratus deck 
to thicken as the inversion and LCL heights diverge. This 
case is obviously quite different from the uncoupled case, 


and it would result in a markedly different forecast. 


C. COUPLED MODEL VARYING MIXED LAYER DEPTH RESULTS 

As noted earlier, this model is useful in that it not 
only provides a 24-hour local forecast but conditions can be 
varied and examined for their effect on the output. A couple 
of "what if" circumstances will be examined for the current 
problem. First, "what if during local maneuvers of the task 
force an ocean front is crossed and the mixed layer is 
suddenly 10 meters deep rather than the current 2 meters?" 
Will this affect the model output? Examining Figure 11, the 
new MLD is evident with the same strong temperature gradient 
as in the previous case. Changes in surface temperature no 
longer have the strong effect on MLD previously noted. This 
1S correct as the heat capacity of a 10 meter mixed layer is 
much greater than that of a 2 meter mixed layer. The near 
surface heating which took place in the previous case late in 
the period is also repeated in this case. A time variance of 
the surface MLD is shown graphically in the lower panel of 
Figure 12. The same wind speed profile as used in the two 
previous cases is evident. Of interest is the large variance 


in the specific humidity profile. 


4l 





ee LURE E TEMPERATURE 


22 18 18 20 22 
2E 
t~ 
a 
| ir م‎ 
2300 4C0 
E luus TEMPERATURE 
18 29 22 
Ir A 
6 
Li 
CH 
1400 1900 
Figure ll. OBL Forecast With 10151351 MLD Set at 1U Meters. 


All Other Input Values are Held Constant 
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Eure 12. Coupled Air Sea Boundary Layer 24-Hour Forecast 
Varying Initial Sea Surface MLD to 10 Meters 
Vice Initial Input of 2 Meters 
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This case looks identical to the uncoupled case with no 
humidity decrease as in the previous case. The air 
temperature profile also matches the uncoupled case while the 
SST has the same trends as in the coupled case. The 
amplitude of variance is much less in this case, and the SST 
and air temperature are never equal. As discussed before, 
heat and moisture are being transferred into the atmosphere. 
However, prior to the two temperatures becoming equal, the 
effects of solar radiation upon the SST cause the two 
temperatures to diverge again. The LCL and inversion plots 
also have the same general characteristics of the uncoupled 
model. The only difference is the slight divergence of the 
two heights late in the period which results in no clouds 


being formed during the period. 


D. COUPLED MODEL VARYING WIND SPEED RESULTS 

The second "what if" case to be examined is one in which 
the wind speed is varied. What if the winds increased from 
Busscurrent conditions to l0 knots late in the period? As 
shown in Figure 13, the extra mixing reduces the SST rapidly 
and drives the MLD down much more rapidly than in the pre- 
vious cases. The early temperature reduction is much stronger 
than in previous cases. There is no near surface heating 
mie ۳ the period which occurred in both of the coupled 


cases examined previously. The rapid decrease in MLD is again 
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Figure 13. OBL 24-Hour Forecast With Wind Speed Increasing 
from 4 to 10 Knots in the Forecast Period 
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shown graphically in the lower panel of Figure 14. : The slow 
increase in wind speed is depicted in the next panel. 

The specific humidity has the same general trends shown 
in Figure 10; however, the gradients are much steeper early 
in the period and fall off rapidly when clouds begin to form. 
The increased wind speed allows the heat and moisture to be 
transferred into the atmosphere much faster than in the 
previous cases. This increased mixing is also apparent in 
the rapid convergence of the sea surface and air tempera- 
tures. Also, the SST shows a continual decrease throughout 
the period as heat is transferred out of the water creating 
convective turbulence in the upper ocean. However, the early 
mation of clouds effectively reduces the incoming short 
wave radiation which would heat the sea surface. 

Cloud top cooling affects the MABL temperature between 6 
and 10 hours into the forecast period. However, this effect 
is negated by the trapping of heat in the boundary layer 
between the stratus deck and the sea surface. This effect is 
apparent during the latter half of the model run. The inver- 
sion height moves above the LCL almost immediately after the 
first increase in wind speed. As the wind increases, the 
depth of the stratus layer also increases. While the initial 
ENOUd top cooling does lower the LCL slightly, the surface 
heating quickly overcomes this effect. This causes the LCL 


to rise late in the period and affects cloud base height. 
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V. RESULTS AND CONCLUSIONS 


While no verification data exist to indicate which of the 
above forecasts was most correct, the primary goal of showing 
the utilization of the model has been demonstrated. As with 
any computer generated product, the forecast generated is a 
rect reflection of the quality of the initialization data. 
The finest computer model will generate poor output given 
poor input. The requirement for more man-machine cooperation 
with this model is such that, with the use of a little common 
sense and some meteorological theory, the program should 
prove useful to the naval geophysicist. Operating in data 
Space regions and often adverse communications areas the 
Meaty to use local conditions as inputs to a locally 
generated forecast should improve forecaster performance. 

In the numerous runs which have been completed the 
performance of the model has proven to be at least a good 
predictor of trends. While often little difference exists 
between the coupled and uncoupled model outputs it is those 
cases which are critical to naval operations that the 
difference is appreciable. Regions of fog and stratus 
formation is one of these circumstances. The formation of 
9 CaN Dé Critical to the ability of naval aircraft being 


A leto accomplish their mission. 





EM/EO propogation is strongly affected by changes in the 
temperature and/or humidity profiles. The ability for a task 
force screen to properly guard a carrier or for a task force 
to remain hidden from enemy radar lies in its abjlityzto 
properly use the environment. Changes in the MLD and the 
subsequent focusing or ducting of sound can be used to find 
enemy targets as well as to hide convoy noise from these same 
forces. 

Improvements in model input techniques and coupling of 
output from this model to IREPS would provide an improved 
package. Making inputs as straight forward and non- 
subjective as possible will aid the fleet operator in 
obtaining a useful product for presentation purposes. Having 
an onboard capability to produce short range single station 
forecasts should help the environmentalists in better serving 
fleet operations. Through proper use of the Air-Sea Boundary 
Layer Model and other environmental data, the trust in 
forecasts presented should improve, and the readiness of 
other fleet units will improve by the efforts of the entire 


geophysics community. 
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